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We discuss phenomenological tests for the frozen infrared behavior of the running coupling 
constant and gluon propagators found in some solutions of Schwinger-Dyson equations of the gluonic 
sector of QCD. We verify that several observables can be used in order to select the different 
expressions of Qs found in the literature. We test the effect of the nonperturbative coupling in the 
r-lepton decay rate into nonstrange hadrons, in the p vector meson helicity density matrix that are 
produced in the Xc2 —> pp decay, in the photon to pion transition form factor, and compute the cross 
sections for eleistic proton-proton scattering and exclusive p production in deep inelastic scattering. 
These quantities depend on the infrared behavior of the coupling constant at different levels, we 
discuss the reasons for this dependence and argue that the existent and future data can be used to 
test the approximations performed to solve the Schwinger-Dyson equations and they already seems 
to select one specific infrared behavior of the coupling. 
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I. INTRODUCTION 

The interface between perturbative and nonperturbative QCD has been studied for many years. To know how far 
we can go with perturbative calculations and how we can match these ones with nonperturbative quantities, obtained 
through sophisticated theoretical or phenomenological approaches, will probably still require several years of confront 
between these methods. However there are phenomenological indications that this connection may not be abrupt 
and the transition from short distance (where quark and gluons are the effective degrees of freedom) to long distance 
(hadron) physics is indeed a smooth one [1]. Actually it has been suggested that the strong coupling constant freezes 
at a finite and moderate value [2] , and this behavior could be the reason for the claimed soft transition. The freezing 
of the QCD running coupling at low energy scales could allow to capture at an inclusive level the nonperturbative 
QCD effects in a reliable way [1,3]. 

The aspects described above reflect a perspective to attack the short and long-distance QCD interface from the 
perturbative side. On the other direction we have nonperturbative methods to investigate the infrared QCD behavior. 
One of these methods, based on the study of gauge-invariant Schwinger-Dyson equations, concluded that an infrared 
finite coupling constant could be obtained from first principles [4] . A list of nonperturbative attempts to determine 
the infrared behavior of the gluon coupling constant and propagator can be found in Ref . [5] , where it is also pointed 
out that the nonperturbative results for these quantities may differ among themselves due to the intricacies of the 
nonperturbative methods and to not always controllable approximations. 

Recently new theoretical results about the infrared behavior of the gluon propagator and the running coupling 
constant appeared in the literature. We have a renormalization group analysis implying in an infrared finite coupling 
[6], and also new numerical lattice studies simulating the gluon propagator, all of them consistent with an infrared 
finite behavior [7]. But the most interesting for us are the new solutions for the gluon and ghost Schwinger-Dyson 
equations that have been obtained with better approximations [8-11]. The absence of certain angular approximations 
in the integrals of the Schwinger-Dyson equations (SDE) are some of the improvements in the Refs. [9,10], which led 
to a value for the infrared coupling constant at the origin (as(0)) roughly a factor 0(3) below the previous ones [12]. 

Assuming that the coupling constant and the gluon propagator are infrared finite we can divide the Schwinger- 
Dyson solutions for the gluon propagator in two classes. One where the gluon propagator is identical to zero at the 
momentum origin [8-11] and another where the propagator is roughly of order l/m^^ [4], where nig is the dynamical 
gluon mass. In Ref. [5] the phenomenological calculation of the asymptotic pion form factor was used to show that 
the experimental data seems to prefer the Cornwall's solution for the coupling and gluon propagator [4]. 

In this work we will extend our approach [5] to others observable quantities, testing the recent nonperturbative 
results obtained for as and for the gluon propagators. These tests span from a purely perturbative calculation to cross 
sections containing some model dependence. We verify that these quantities are affected by the infrared behavior 
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of the coupling in different amounts. However we can certainly see that the existent data can be used to test and 
distinguish different nonperturbative calculations of the infrared QCD behavior. 

There arc several important points to stress. The SDE solutions represent attempts of a very ambitious program 
aiming to obtain the behavior of the coupling constant at all scales. Of course, these equations contain approximations 
and there are even questionings about the definition of the coupling in the infrared, but once these facts are assumed 
and knowing that in Nature there must be only one coupling constant we should be able to compare these results 
to the data. Some of the tests we propose here are purely perturbative, but they still can say something about 
the approximations used to obtain as- Another point is that the whole idea of a running charge is that it is to be 
used in a skeleton expansion of graphs representing a process, an expansion which automatically subsumes infinitely 
many graphs, in a scheme emphasized by Brodsky and collaborators [13] . But if, contrary to many phenomenological 
evidences [1,3,13], the infrared value of the coupling constant were large enough (how large depends on the process), 
it is hard to sec how one can make sense out of inserting such a large value into perturbative expressions. This point is 
at the heart of the discussions that we shall present in this paper. The SDE solutions are determined forcing its high 
energy behavior to be identical to the perturbative one, after this is done the infrared behavior is totally dictated by 
the equations. This procedure is a natural one because at large momenta the coupling constant is perturbative and 
the QCD determination of the leading S matrix terms is meaningful, although the perturbative quantities depend on 
an infrared paramotrization hidden in the QCD scale ( Aqcd)- 

The solutions that we shall discuss have different values for as{0). It could be argued that the value of Qs(0) for one 
of these SDE solutions is too small to explain chiral symmetry breaking in the usual model with essentially one-gluon 
exchange [14]. It has been known for many years that this sort of model of chiral symmetry breaking doesn't work, 
precisely because it does require such large values for 0^(0). In QCD, unlike QED, chiral symmetry breaking is easily 
explained -indeed, required by- confinement, as numerous authors in the eighties have explained in detail. A value 
of as{0) of order one or larger would be grotesquely out of line with many other phenomena, such as non-relativistic 
potentials for charmonium [15]. Finally, this is also what comes out when performing an analytic perturbation theory 
[16]. Moreover, the fact is that the nonperturbative solutions do not have only different as{0) values but they also run 
differently near the infrared. The difi^erent infrared parametrization and running will leave a signal in the comparison 
with the experimental data. 

The distribution of the paper is the following: In Section II we present a comparison between the most recent 
calculations of the infrared behavior of the QCD running coupling constant. Section III compares the nonperturbative 
value of as with its value measured in r-lepton decay into nonstrange hadrons. In Section IV we discuss the eff'ects of 
a frozen coupling constant in the diagonal elements of the p vector meson helicity density matrix that are produced 
in the Xc2 PP decay. Section V is devoted to a discussion of the 7 ^ ttq transition form factor. In Section VI we 
revisit the calculation of Ref. [17], where the Pomeron model of Landshoff and Nachtmann for elastic proton-proton 
scattering was used to restrict the value of the dynamical gluon mass, in the light of the new SDE solutions. The 
same model is used in Section VII to compute the exclusive p production in deep inelastic scattering. Section VIII 
contains our conclusions. 



II. INFRARED BEHAVIOR OF THE RUNNING COUPLING CONSTANT 



In this Section we present the results for the running coupling constant obtained recently through the solution of 
SDE for the gluon and ghost sectors of QCD. As remarked in Ref. [5] and in our introduction, it is important to stress 
that the different solutions appear due to the different approximations made to solve the SDE, which, unfortunately, 
are necessary due to their complicated structure. 

The most recent calculations of the as infrared behavior do not make use of an approximate angular integration 
in the SDE, resulting in a value of as{0) a factor of 0(3) below the former results. This is the case of the as{q'^) 
determined by Fischer and Alkofer which is given by [9] 

/X Q:a(0) 

"^^^''^~ln(e + aia;«2 + 6iar''2)' 



where 



0^(0) = 2.972, 

ai = 5.292 GeV-^"^ 

02 = 2.324, 

61 = 0.034 GeV-^''% 

62 = 3.169. 
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We also have an ansatz for the running coupHng proposed by Bloch. Its expression for as{q^) has the following fit 



[8]: 



asB{q^) = a{lK^QCD) 
1 



(2) 



+ 



coao + 



47r 



/?o Vlog(0 /-I 



where I = q /Aqqj^, cq = 15, ao = 2.6, and /3o = 11 — 3"/, where n/ is the number of flavors. 

Equation (1) has been obtained solving the coupled SDE for the gluon and ghost propagators in the Landau gauge. 
Eq. (2) is an ansatz motivated by earlier results done with the help of the angular approximation. Both are consistent 
with a gluon propagator that vanishes in the infrared. The difference in the couplings, which is clearly displayed in 
Fig.(l), is that Bloch fix the scale of the coupling constant at Kqcd = SOOMcV while Alkofcr and Fischer [9] fix it 
comparing to the experimental input at as{Mz), unfortunately their calculation when compared to a perturbative 
expression for coupling constant in the MOM scheme lead to a value Aqcd = 715MeV which seems to be too large. 
The difference between the values of as{0) is negligible, while the difference at intermediate momenta is going to be 
transferred to the phenomenological quantities that we will compute. 




FIG. 1. Comparison between the running couplings obtained from different solutions of Schwinger-Dyson equations. 



The other expression for the running coupling constant that we will use in the next sections is the one determined 
by Cornwall many years ago [4] 



asciq'^) 



- 



/3oln [{q^+AM^{q^))/A^] 



(3) 



where Mg{q^) is a dynamical gluon mass given by, 



M^{q^) = ml 



ln(^^^^ 
ln(%i) 



-12/11 



(4) 



A(= Aqcd) is the QCD scale parameter. This solution is the only one that has been obtained in a gauge invariant 
procedure. The gluon mass scale has to be found phenomenologically. A typical value is [4,17] 
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rUg = 500 ± 200 MeV 



(5) 



for A = 300 MeV. The Bloch and Cornwall's expression at large momenta map perfectly into the perturbative running 
coupling, as can be seem in Fig.(l). We will refer to these different results as model A, B and C respectively. 



III. T-LEPTON DECAY RATE INTO NONSTRANGE HADRONS 



The tests of the running coupling constant behavior that we shall discuss will depend on the exchange of gluons at 
different levels, i.e. the softer is the exchanged gluon the more we can test the nonpertubative behavior of a^. The 
physics of the r-lepton decay rate into nonstrange hadrons is one where QCD can be confronted with experiment to a 
very high precision. The measurement of ag in these decays and its comparison with the nonperturbative expressions 
will correspond to what can be called our most perturbative test. 

The normalized r-lepton decay rate into nonstrange hadrons (/is=o) is given by [18] 

r(T ^ hs=oiy) 

= -^c|Kd| Sew{'i- + 6ew + Sp + 6Np), (6) 

where A^c is the number of colors, the flavor mixing matrix element is \Vudf = 0.9511 ± 0.0014 and Sew = 1.0194 
and Sew = 0.001 are electroweak corrections. The nonperturbative corrections are small and consistent with zero, 
Snp = —0.003 ± 0.003. In particular, these corrections are not directly related to the nonperturbative infrared 
behavior of as- The term Sp represents the perturbative QCD effects and its value can be calculated if we have 
RiT)\s=o- With the experimental result obtained by the ALEPH Collaboration [19], i?(T)|^!fo*=3.492±0.016, and the 
other values above one can estimate the perturbative corrections: 

6p^* = 0.203 ± 0.007. (7) 

On the other hand dp can be calculated in the framework of perturbative QCD and, in the MS scheme, is given by 
the expansion 

J*? = + 5.2023 + 26.366 (^)' + • • • , (8) 

where is the running coupling constant taken at the scale of the r-lepton mass M^. = 1.777GeV. Substituting 
Eq.(7) in the left hand side of Eq.(8) and solving the resulting polynomial equation in one obtain 

^expt ^ Q 34Q4 ^ 0.0073. (9) 

Now we can compare the above experimental result with the different values of the nonperturbative running coupling 

at the r-lepton mass scale. Some comments arc in order before such comparison is performed. The nonperturbative 
couplings have been obtained in the case of zero flavors. When we compare the perturbative to the nonperturbative 
results it should be understood that the introduction of fermions will affect the three models in the same direction 
(increasing as(0)) and this effect would be small for a small number of flavors [4,5]. Another point is that Eq.(8) 
has a renormalization scheme dependence, while the nonperturbative solutions have been obtained in one unrelated 
approach. However the nonperturbative couplings of the three models are obtained forcing the matching with the 
perturbative running coupling at very large momenta and at leading order there is no problem with scheme dependence, 
in such a way that the comparison can be made but not with the precision of the experimental result. Using Eqs. (1) 
to (3) (with rUg = 500 MeV and A = 300 MeV) instead of the perturbative running coupling we obtain 

asA = 0.676 , asB = 0.354 , Usc = 0.348. (10) 

As one could expect only the last two results are compatible with Eq.(9), and the model A fail at a level much above 
the one where we could claim that any scheme dependence is important. This shows the importance of the scale 
used to fix the coupling constant. Note that Alkofer and Fischer [9] (model A) constrain the coupling comparing its 
value to the experimental one of as{Mz), which seems to be the most rigorous way of doing it, unfortunately their 
calculation lead to a value Aqcd = 715MeV, which is too high when compared to most of the values found in the 
literature. As the r-lepton mass scale is one where we believe that perturbative QCD can already be applied we have 
to conclude that the SDE solutions of model A does not describe the experimental data at such energies. Model B 
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docs it quite well but it was fitted to a smaller Aqcd value. Note that the model A docs not fit the data not because 
it has a large value for as{0), which is not so difl:erent from the model B, but because the running coupling of model 
A decreases much more slowly. 

As we claimed previously, this test is the most perturbative one in the sense that it does not involve the introduction 
of hadronic wave functions directly in the calculation. In the next sections we always will need to introduce some extra 
hadronic information besides the coupling and giuon propagator. It is important to stress that this result does not 
settle the question about which one is the correct solution of gluonic SDE, it only give hints to which approximations 
may be or not suitable to be performed when solving these equations. 



IV. EFFECTS OF A FROZEN COUPLING CONSTANT IN THE AMPLITUDES OF THE xc2 ^ pp 

DECAY 

It is known, for a long time, that a dynamical gluon mass can affect quarkonium decays [4,20,21]. In the study of 
Refs. [20,21] it can be noted that the dynamical mass effect is relevant in the study of their hadronic decays mostly 
due to phase space factors. Therefore not all quarkonium decays are good candidates to probe the infrared behavior 
of the running coupling constant, because the effects of phase space overwhelms the changes in the running coupling. 

It is also known that the exclusive quarkonium production in pp or pp interactions and subsequent decay is strongly 
dependent on the internal structure of the hadrons involved, which are described by their hadronic wave functions. In 
the inclusive case the important dependence is on the distribution function. Perturbative QCD fixes the asymptotic 
form of these functions as ^ oo and their general evolution, but in any realistic calculation they have to be taken 
as phenomenological quantities, to be experimentally determined via a set of physical information and then used in 
other processes. 

Even when we are dealing with heavy quarks, and try to propose tests for the interface between perturbative 
and non-perturbative QCD, we verify that we cannot test the infrared behavior of quarks and gluons vertex and 
propagators independently of hadronic wave or distribution functions. In this Section we discuss the effects of a 
frozen coupling constant in the measurement of the p vector meson helicity density matrix that are produced in the 
Xc2 PP decay. This is a case where the infrared behavior of the coupling is entangled with the behavior of the wave 
or distribution functions. We advance that in Section V we will discuss a more ordinary case where the wave function 
is well known, and where it is possible to study the QCD infrared behavior with less uncertainty [17]. 

Some time ago Anselmino and Murgia [22] considered the Xc2 pp decay process of polarized charmonium states 
created in pp or pp interactions and have shown how the observation of the polarization of the vector meson, via a 
measurement of its diagonal helicity density matrix elements, neatly depends on the p wave function and helps in 
discriminating between different kinds of these quantities. 

The study of Ref. [22] provide an interesting arena to introduce the effects of a frozen running coupling constant. 
Our result will show that depending on the form of the wave function we could observe a larger or smaller effect of 
the freezing in the infrared, showing how the nonperturbative behaviors of wave (or distribution) functions with the 
ones of the effective coupling can become entangled. 

The processes that we consider are the exclusive 

PP^Xc2^PP (11) 

or the inclusive 

pp^Xc2 + X ^ pp + X (12) 
production of a pair of p vectors with the subsequent decay 

p ^ TTTT, (13) 

and the quantity experimentally observed is the angular distribution of either one of the pions in the helicity rest 

frame of the decaying p. 

The pion angular distribution depends of the spin state of the p via the elements of its helicity density matrix 
P\,X'{py- 

3 

W(e,$) = — [po.ocos^ 9 + (/9i,i - pi-i) sin^Gcos^ $ + (pi,i + pi,_i) sin^ 6 sin^ $ 

-\/2(iiepi,o)sin2ecos$], (14) 



5 



where O and <[> are, respectively, the polar and azimuthal angles of the pion as it emerges from the decay of the p, 
in the p helicity rest frame. Eq. (14) can be integrated over or $ generating polar and azimuthal distributions, 
which can be measured and give information on px,\' (p)- The details of this procedure (i.e. the vector meson helicity 
density matrix calculation for massless quarks) can be found in Ref. [22] and references therein. 
The helicity density matrix of the p meson are [22] 

px,x'=0, A^A', po,o = 1 - 2/91,1 (15) 



= = ^ r , (16) 

^1 + 3 LA" "i^ F(e) 

|Ai,-iP 

where the reduced amplitudes Aa.a' are given in Eqs. (2.10), (2.11) and (2.14), (2.15) of Ref. [23] and do not depend 
on the p production angles 9 and 0, but do depend on the wave functions. Eq. (16) has the same form both for 
exclusive and inclusive p production, but the dependence on the production angle 9 is different in the two cases 

pexfn-. _ COS^(^) ^ 7^ 



^ ' 1 + 6COS2(0)+COS4(0)- ^ ' 

The ratio of reduced amplitudes |Ao,oP/|Ai,_ip appearing in Eq.(16) are computed as a function of the p longitudi- 
nally (L, A = 0) and transversely (T, A = ±1) polarized vector mesons wave functions, which are indicated respectively 
hy ipL and <^t, and correspondent decay constants given by /l and /t- This ratio is equal to [22] 

0,01 ^ h] l-'o,o| ^^g^ 



|Ai _i] V6 KhJ \h, 
where 



T ^ [\ ^ ( (\\ ( A^ ^^(^^^xH(i-a:)(i-y)M^] 



Io,o = - ^ 



dxdyipL{x, Q^)(/3l(j/, Qy) 



.2 a,{xyMl)a[{l - x){l - y)M^ 



xy{l -x){l- y){2xy -x-y) 



{x - yf 



2xy -x-y 



(21) 



In the above equations Af^ = 3.6 GeV is the Xc2 mass and Q^. — min(a;, 1 — x)Q. 

To calculate the helicity density matrix pi^i we choose two different sets oi^fL and ^pT'- a set of symmetric distribution 
amplitudes 



and the QCD sum rules amplitudes [24] 



^L = ^ 6x'(l - x) 



ipL{x, Q;) = 6a:(l -x){l+ ^-Cl'\2x - 1) 



~ 2 " 

a,(Q^) 


2/3 1 







(22) 



(23) 



a2 



frix^x) = 6a;(l - x) 



"^(Mt) 



4/25 



lcl'\2x-l) 



_2 -| 52/75' 



(24) 



where /Lt|, = 0.5 (GeV/c)^, /xf. = 0.25 (GeV/c)^, and denotes Gegembauer polynomials. In both cases we assume 
fL = fx- 
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FIG. 2. Values of the helicity matrix element as a function of the p-meson production angle 9. The exclusive process is 
computed with the symmetric distribution amplitudes for the Bloch and Cornwall's coupling constants. 



1 — I — 1 — I — 1 — I — 1 — I — 1 — I — ' I ' I ' I ' I ' I 




1 — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I 
20 40 60 80 100 120 140 160 180 200 

e(deg) 



FIG. 3. Same notation as in Fig. (2) but calculated with QCD sum rules distribution amplitudes 

Our results for the exclusive process (Eq.(ll)) are displayed in Fig.(2) and Fig. (3). Fig. (2) shows pi^i as a function 
of the /5-meson production angle 9, computed with the symmetric wave function whereas Fig. (3) was calculated with 
the QCD sum rule amplitudes. Both figures show curves for different behaviors of the running coupling constant, the 
one of Eq. (2) obtained by Bloch and the one of Eq. (3) obtained by Cornwall with nig ~ 600 McV and A = 300 McV. 
The calculations for model A are not shown in the figures. The results for this model are similar to the ones of model 
B with curves located in the opposite direction to the curves of model C. 
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FIG. 4. Values of the holicity matrix clement pi,i as a function of the p- meson production angle 6 in the case of inclusive Xc2 
production. The curves were obtained using the symmetric distribution amplitudes and the expressions of Eq.(2) and Eq.(3) 
for the coupling constant. 
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FIG. 5. Same notation as in Fig. (4) but calculated with QCD sum rules distribution amplitudes 

As claimed by Anselmino and Murgia [22] the measurement of the hehcity matrix element pi^i of the vector meson 
in the Xc2 PP can indeed discriminate differences between wave (or distribution) functions, but the effect of different 
running coupUngs is not negligible and one effect can mask the other. Of course, this is not an easy experiment, but 
is a feasible one. If the measurement is performed with high precision it can provide information about the infrared 
behavior of the coupling constant. 

The curves for the inclusive process (Eq.(12)) are displayed in Fig. (4) and Fig. (5). Fig. (4) shows pi^i as a function of 
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the p- meson production angle 9, computed with the symmetric distribution amphtude whereas Fig. (5) was calculated 
with the QCD sum rule amplitudes. The differences in the curves are of the same order as in the exclusive case. 

Note that in both cases (exclusive or inclusive) the values of pi^i is modified by a factor of 0(2 — 3) when we use 
the wave functions obtained through QCD sum rules. This is expected because they have a stronger dependence 
on as{q'^). The difference occurs at different angles for the exclusive and inclusive Xc2 production. Actually, it is 
surprising that in such a delicate experiment, where the effects of the running coupling appear in a ratio like the one 
described by Eq.(19), there is still a clear signal of the dependence on the infrared value of the coupling constant. 
In this section we verified how the infrared behavior of the coupling can be masked by the nonperturbative wave 
functions. In the next sections we will look at processes where the effects of wave or distribution functions are not so 
pronounced or are better known. 



V. THE 7 ^ TTo TRANSITION FORM FACTOR 



When making predictions within pcrturbativc QCD, we are always confronted with problems like the choice of 
distribution amplitudes as discussed in the previous section, the choice of renormalization scale /x and scheme for 
the coupling constant, etc... However there are calculations that have been discussed extensively in the literature, 
processes where the measurement is more sensitive only to the asymptotic distribution amplitude of a particular 
hadron, and where an optimal renormalization scale has been estimated. This was the case of the pion form factor 
discussed in Ref. [5], and is the case of the 7 ^ ttq transition form factor to be described here. 

The photon-to-pion transition form factor Fj„{Q'^) is measured in single-tagged two-photon e+e" e+e~7r° 
reactions. The amplitude for this process has the factorized form 

4 

F^^{Q^) = -^j^ dx4>^{x,Q^)T^,{x,Q^), (25) 
where the hard scattering amplitude T^^{x,Q^) is given by 

Using an asymptotic form for the pion distribution amplitude = V3fnx{l — x) we obtain [25] 

Q%AQ') = V.{i-l^) (27) 

where Q* = exp~^/^ Q is the estimated Brodsky-Lepage-Mackenzie scale for the pion form factor in the scheme 
discussed in Ref. [25] . 

In Fig. (6) we compare the photon to pion transition form factor with CLEO data [26]. The curves were computed 
with different expressions for the infrared behavior of the running coupling constant. We assumed /^r ~ 93 MeV 
and A — 300 MeV. Using the running coupling constant given by the expression of Eq.(2) we obtain a fit for the 
photon-pion transition form factor very far from the experimental data. The result obtained when we use Eq.(l) 
is not shown and gives an even worse fit. The infrared value of the coupling constant is so large in the case of the 
coupling constants given by Eqs.(2) and (1), that we are not sure that the perturbative result can be trusted even at 
such large momentum scale. 
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FIG. 6. The 7 — > tt transition form factor calculated with different expressions for the infrared behavior of the running 
coupling constant. 

The infrared coupling eonstants related to the class of SDE solutions consistent with an infrared finite propagator 
that vanishes at origin of momenta, are much stronger than most of the phenomenological estimates of the frozen 
as{0) value that we quoted in Ref. [5] (as(0) « 0.7 ± 0.3), and are at the origin of the strange lower curve of Fig. (6). 
Notice that due to the large scale of momenta involved in this process it is not expected that higher order corrections 
are still important. The data is only compatible with Eq.(3), which has a smoother increase towards the infrared 
region. Perhaps this behavior is actually indicating that the transition to the infrared should be a soft one. Again 
the result for model A is not shown because it is even worse than the one of model B. 

Note that in Fig. (6) the curves obtained with the Cornwall's coupling constant do not show large variation in the full 
range of uncertainty of the dynamical gluon mass (see Eq.(5)). It is interesting that its behavior is quite stable in this 
case as well as for the pion form factor studied in Ref. [5] . If we had large variations of the infrared coupling constant 
with the gluon mass scale we could hardly propose any reliable phenomenological test for its freezing value. We stress 
that the results are obtained for a perturbative scale of momenta and it seems that we have to choose between two 
possibilities: This process cannot be predicted by perturbative QCD up to a sc;ale of several GeV, or some of the ESD 
solutions are predicting a too large value of the coupling constant in the infrared and the approximations made to 
determine these solutions are too crude. 



VI. ELASTIC DIFFERENTIAL CROSS SECTION FOR PP SCATTERING 



We have verified in Section IV one case where the behavior of the coupling constant in the infrared is entangled with 
the one of the distribution amplitudes. In Section V it was shown a classic QCD calculation where the asymptotic 
behavior of the wave function is more important than the infrared one, leading to a cleaner test for the nonpertur- 
bative behavior of the coupling constant. Here we will discuss a test for the infrared QCD behavior that makes use 
of a phenomenological model for diffractive interactions [17]. In this case we will also make use of a hadronic wave 
function, but the main source of uncertainty is in the fact that we assume a model for elastic scattering containing a 
series of approximations, which are beyond the scope of a pure QCD calculation. 

We will compute the elastic differential cross section for pp scattering within the LaiicLshoff-Nachtmann model 
(LN) for the Pomeron exchange [27]. As discussed in Refs. [17,28], this is a model where the Pomeron is represented 
by two-gluons exchange, and it is particularly dependent on the infrared properties of the gluon, not only on the 
coupling constant but also on the behavior of the gluon propagator. In this model one of the gluons carry most of 
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the momentum exchanged in the interaction while the other seems to appear just to form the colorless Pomeron. 
Therefore besides the coupling constant we will also need the expressions for the gluon propagators obtained through 
the solutions of the SDE, although this dependence with the nonperturbative behavior of the gluon propagator is 
important only when the gluon is exchanged at low momentum. The gluon propagator in Landau gauge is written as 



where the expression for D{q'^) obtained by Cornwall is given by 



q'^+iM^ 
A2 



(28) 



(29) 



We will also make use of the running coupling constant obtained by Fischer and Alkofer (model A) and their respective 
propagator D{q^) = Z{q'^)/q'^, where Z{q^), in Landau gauge, is fitted by 



Z{x) 



I W l+2<5 



and 



R{x) = - 



dx^'^ 



(30) 



(31) 



where the constants appearing in Eq.(30) and Eq.(31) are given by 
aAin"^) = 0.9676, 
K = 0.5953, 
6 = -9/44, 
c= 1.8934 GeV"^'', 
d = 4.6944 GeV"^'". 
In the LN model the elastic differential cross can be obtained from 

da _ \Ais,t)\'^ 
H ~ 167rs2 

where the amplitude for elastic proton-proton scattering via two-gluon exchange can be written as 

A{s,t)=isB>al[Ti-T2] 

with 

Ti = y d^kD (I + fc) Z) (I - k) \G^{q, 0)|2 

T2 = j d^kD (I + fc) £) (I - k) Gp (q, fc - I) [2Gp(g, 0) - (q, k - |)" 
where Gp{q, k) is a convolution of proton wave functions 

GM,k) = jd^vdnri^,p)i,i.,V-k-.q). 

Gp{q, 0) is given by the Dirac form factor of the proton 

/ ^ ^ / „x 4to2 - 2.79t 1 
F.it) = Gpiq,0)= (i_,/o.71)2 - 

To estimate Gp{q, k — q/2) we assume a proton wave function peaked at k = 1/3 and obtain 



Gp(q,k-^)=F,(^q^ + 9 



(32) 



(33) 



(34) 



(35) 



(36) 



(37) 



(38) 
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Eq.(32) to Eq.(38) can be computed using the couplings and propagators discussed above. We compare the differen- 
tial elastic cross-section for proton-proton scattering with the experimental data of Breakstone et al. at y/s = 53 GeV 
[29]. The results have to bo adjusted by a normalization factor gO i^* which accounts for the energy dependent part 
of da/dt [30]. For large \t\ values, double- Pomeron exchange and three-gluon exchange are known to be important; 
we therefore do not expect to describe the data near and above -t ~ 0.5 GeV"^ [30,31] 

The data of elastic proton-proton scattering is well fitted assuming a dynamical gluon mass of nig = 370 MeV 
for A = 300 MeV. The calculation has a small variation with the value of TUg and is more sensitive to the ratio 
nig/ A. Using the coupling and propagator of Alkofer et al. (model A) we obtain a curve that is about one order of 
magnitude away from the experimental points. Unfortunately we do not have the propagator solution for model B 
and for this reason we do not show the results in this case, although we can guess that it would produce a curve in 
Fig. (7) between the curves of model A and C. 




0.01 ' I 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 

-t(GeV') 

FIG. 7. Differential pp elastic cross section at y^s) = 53 GeV computed within the Landshoff-Nachtmann model for the 
Pomeron, using different infrared couplings and gluon propagators obtained from SDE solutions. 

There is a striking difference between the two classes of SDE solutions that we discussed in the Introduction. In 
the case of Cornwall's solution it is the product g^D{q^) that behaves roughly as '^/mg^ as ^ and this product 
has no dependence [4] . This docs not happen for the class of gluon propagators that seems to vanish at origin. As 
in this case the coupling is not canceled in the da/dt calculation and as it is a factor of 0(3) larger, it is not difficult 
to understand the one order of magnitude difference in the result of the cross sections. This specific calculation is 
model dependent, however it is quite sensitive on the infrared behavior of the theory. 



VII. EXCLUSIVE p PRODUCTION IN DEEP INELASTIC SCATTERING 

Donnachie and Landshoff [32] successfully described the process "{*{q)p{P) p{r) p{P') with a soft Pomeron 
exchange, and later they [33] considered Pomeron exchange as being the exchange of two gluons as in the previous 
section. Pursuing that idea Cudell [34] proposed the following expression for the polarized differential cross section 
for exclusive p production: 

^ = l^.l'*') [3^iW]' , j=T, L, (39) 

where agim — 1/137 is the electromagnetic coupling constant and w"^ = {q + P)^. The factor $ is given by 
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where fp ~ 30MeV is the p form factor and mp ~ 770 MeV is the p mass. 

/ 2 \ 0.08+a't 

Wo 

with Wq = 1/a' ~ 4GeV^. The Dirac form factor of the proton is 

Ami - 2.79t 1 



Fiit) = 



Aml-t (l-i/0.71)2' 



where mp is the proton mass. 
The ampUtude is 



(41) 



(42) 



p,ll=ml + Q\ j = T,L 

where Pt = Pj^ = Prp ■ (m^ + + t)/{2mpQ). In the last term of the integrand we have the frozen couphng 

constant a„ = as{0) and the gluon propagator D{p'^). 



10'- 



10° 



10=- 



J3 



D 10 

loH 

10'- 



^ 

Total cross section exclusive p production 
-■-CORNWALL mg=370 MeV 
-•-ALKOFER etal 
• EMC & LAME experimental data 
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Q'[GeV'] 



FIG. 8. Total cross section for exclusive p production in deep inelastic scattering. 
The total cross section can be obtained through the integral 

atotaiiQ') =ldt{^+ ' ^ « (43) 

This cross section is dependent not only on the nonpcrturbativc coupling constant but also on the gluon propagator 
as the calculation of the previous section. We compared the solutions of Alkofer et al. (model A) and Cornwall (model 
C). The reason for not using model B is the same as the one discussed in the previous section. With the running 
coupling constants (cqs. (1) and (3)) and their respective propagators (cqs. (30) and (29)) we obtain the results shown 
in Fig. (8), where they are compared with the experimental data from EMC [35] and LAME [36]. As explained in Ref. 
[34] the model for exclusive p production is consistent with experiment only in the range of momenta shown in the 
figure, at larger values of there are other contributions to the cross section. 



13 



VIII. CONCLUSIONS 



In this work we proposed phenomenological tests for the infrared behavior of running coupUng constant that have 
been obtained through the nonpcrturbativc sohitions of Schwingcr-Dyson equations of the QCD gluonic sector. We 
are only considering SDE solutions that are infrared finite. Although infrared divergent solutions have not been 
fully discarded, there are several indications that these QCD Green functions are well behaved in the infrared. 
Unfortunately due to the complexity of the SDE, approximations are still necessary to obtain these solutions and 
they lead to different expressions. Phenomenological tests are important because they are able to indicate which are 
the correct approximations that can be performed to solve SDE, selecting the solution that is compatible with the 
experimental data. 

We studied the effects of a frozen coupling constant in the r-lepton decay rate into nonstrange hadrons. This 
is a typically perturbative test and, in particular, it shows the importance of the method used to fix the scale of 
the running coupling constant, which is different for the couplings shown in Eqs.(l) and (2) (model A and B). The 
measurement of the p vector meson helicity density matrix that are produced in the Xc2 — > PP decay is one interesting 
experiment to be performed. The measurement of the diagonal matrix clement provide information on the infrared 
behavior of the coupling constant as well as on the distribution amplitude. In this case both behaviors are entangled 
and a good knowledge on the distribution amplitude is needed in order to select a preferred behavior for the coupling 
constant. The experiment is not an easy one, but can be performed and will give information complementary to the 
others that we have discussed. 

The effects of the infrared behavior of the coupling constant in the photon-to-pion transition form factor F^^{Q'^) is 
much more interesting. The calculation probes the asymptotic behavior of the pion distribution amplitude and provide 
a cleaner test for the freezing of the coupling constant. The differences between the existent solutions obtained from 
SDE for asiq^) are shown clearly in the comparison with the experimental data, which seems to select the coupling 
constant compatible with dynamical generation of a gluon mass. 

The analysis of the differential cross section for proton-proton scattering and the cross section for exclusive p 
production in deep inelastic scattering arc also miambiguously selecting one specific behavior for the QCD Green 
functions in the infrared. In this case we made use of a model of diffractive scattering proposed by Landshoff and 
Nachtmann and it could be argued that the model has not been obtained from straightforward QCD calculations. 
However the model has a great success in the explanation of diffractive physics at low energy, and again the comparison 
with experimental data is consistent with the result obtained from the other phenomenological tests. Diffraction in 
this model is explained through the exchange of two gluons, where one of the gluons carries most of the momentum 
and the other is soft. Therefore the model actually probes the infrared region. 

It is worth mentioning the influence of gauge invariance and presence of fermions in the solutions for the coupling 
constants that we have discussed. The Cornwall's solution was obtained in a gauge invariant procedure, the others 
were obtained in Landau gauge. We do not expect that gauge invariance introduce a large effect in the result. In 
general solutions of SDE are relatively stable in respect to changes in the gauge choice, and we should expect that 
if the solutions minimize the vacuum energy the gauge dependence should disappear. The effect of the number of 
flavors in Cornwall's solution [4] is not so strong, and it appears in the coefficient /?o of the coupling constant and 
in the gluon mass equation increasing the value of the frozen coupling. If a nonzero number of flavors produces any 
observable effect, this one should act in the same sense for all solutions. Therefore, we do not expect large changes in 
our results with the inclusion of fermion loops in the SDE solutions. 

We conclude pointing out that Schwinger-Dyson equations provide a powerful tool to investigate the QCD infrared 
behavior. These are complicated equations whose solutions are determined only after some specific approximations 
leading to different expressions. These expressions can be phenomenologically tested, indicating which are the most 
reliable approximations. In the tests presented here only one solution, consistent with a dynamically generated mass 
for the gluon, was shown to be compatible with the experimental data. 
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